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RESUMEN 

El resumen sera traducido al espanol por los editores. In popular cold dark matter cosmological 
scenarios, stars may have first appeared in significant numbers around a redshift of 10 or so, as the gas within 
protogalactic halos with virial temperatures Tvir ^ 10^"^ K (corresponding to masses comparable to those 
of present-day dwarf ellipticals) cooled rapidly due to atomic processes and fragmented. It is this 'second 
generation' of subgalactic stellar systems, aided perhaps by an early population of accreting black holes in their 
nuclei, which may have generated the ultraviolet radiation and mechanical energy that ended the cosmic "dark 
ages" and reheated and reionized most of the hydrogen in the universe by a redshift of 6. The detailed history 
of the universe during and soon after these crucial formative stages depends on the power-spectrum of density 
fluctuations on small scales and on a complex network of poorly understood 'feedback' mechanisms, and is one 
of the missing link in galaxy formation and evolution studies. The astrophysics of the epoch of first light is 
recorded in the thermal state, ionization degree, and chemical composition of the intergalactic medium, the 
main repository of baryons at high redshift s. 

ABSTRACT 

In popular cold dark matter cosmological scenarios, stars may have first appeared in significant numbers 
around a redshift of 10 or so, as the gas within protogalactic halos with virial temperatures Tvir ^ 10^'^ K 
(corresponding to masses comparable to those of present-day dwarf ellipticals) cooled rapidly due to atomic 
processes and fragmented. It is this 'second generation' of subgalactic stellar systems, aided perhaps by an 
early population of accreting black holes in their nuclei, which may have generated the ultraviolet radiation and 
mechanical energy that ended the cosmic "dark ages" and reheated and reionized most of the hydrogen in the 
universe by a redshift of 6. The detailed history of the universe during and soon after these crucial formative 
stages depends on the power-spectrum of density fluctuations on small scales and on a complex network of 
poorly understood 'feedback' mechanisms, and is one of the missing link in galaxy formation and evolution 
studies. The astrophysics of the epoch of first light is recorded in the thermal state, ionization degree, and 
chemical composition of the intergalactic medium, the main repository of baryons at high redshifts. 
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1. INTRODUCTION 

At epochs corresponding to z ~ 1000 the inter- 
galactic medium (IGM) is expected to recombine and 
remain neutral until sources of radiation and heat de- 
velop that are capable of reionizing it. The detection 
of transmitted flux shortward of the Lya wavelength 
in the spectra of sources at z > 5 implies that the 
hydrogen component of this IGM was ionized at 
even higher redshifts. The increasing thickening of 
the Lya forest recently measured in the spectra of 
SDSS z ~ 6 quasars (Becker et al.2001; Djorgov- 
ski et al.2001) may be the signature of the trailing 
edge of the cosmic reionization epoch. It is clear 
that substantial sources of ultraviolet photons and 
mechanical energy, like young star-forming galaxies, 
were already present back then. The reionization of 
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intergalactic hydrogen at z > 6 is unlikely to have 
been accomplished by quasi-stellar sources: the ob- 
served dearth of luminous optical and radio-selected 
QSOs at z > 3 (Shaver et al.l996; Fan et al.2001), 
together with the detection of substantial Lyman- 
continuum flux in a composite spectrum of Lyman- 
break galaxies at {z) = 3.4 (Steidel et al. 2001), both 
may support the idea that massive stars in galac- 
tic and subgalactic systems - rather than quasars - 
reionized the hydrogen component of the IGM when 
the universe was less than 5% of its current age, and 
dominate the 1 ryd metagalactic flux at all redshifts 
greater than 3. An episode of pregalactic star forma- 
tion may also provide a possible explanation for the 
widespread existence of heavy elements (like carbon, 
oxygen, and silicon) in the IGM. There is mounting 
evidence that the double reionization of helium may 
have occurred later, at a redshift of 3 or so (see Kriss 
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et al. 2001, and references therein): this is hkely due 
to the integrated radiation emitted above 4 ryd by 
QSOs. 

Estabhshing what ended the dark ages and when 
is important for determining the impact of cosmolog- 
ical reionization and reheating on several key cosmo- 
logical issues, from the role reionization plays in al- 
lowing protogalactic objects to cool and make stars, 
to determining the thermal state of baryons at high 
redshifts and the small-scale structure in the tem- 
perature fluctuations of the cosmic microwave back- 
ground. Conversely, probing the reionization epoch 
may provide a means for constraining competing 
models for the formation of cosmic structures: for ex- 
ample, popular modifications of the CDM paradigm 
that attempt to improve over CDM by suppressing 
the primordial power-spectrum on small scales, like 
warm dark matter (WDM), are known to reduce the 
number of collapsed halos at high redshifts and make 
it more difficult to reionize the universe (Barkana 
et al.2001). In this talk I will summarize some re- 
cent theoretical developments in understanding the 
astrophysics of the epoch of first light and the im- 
pact that some of the earliest generations of stars, 
galaxies, and black holes in the universe may have 
had on the IGM. 

2. REIONIZATION BY MASSIVE STARS 

In a cold dark matter (CDM) universe, struc- 
ture formation is a hierarchical process in which 
non linear, massive structures grow via the merger 
of smaller initial units. Large numbers of low- 
mass galaxy halos are expected to form at early 
times in these popular theories, leading to an era 
of reionization, reheating, and chemical enrichment. 
Most models predict that intergalactic hydrogen was 
reionized by an early generation of stars or accreting 
black holes at z = 7 — 15. One should note, how- 
ever, that while numerical N-body+hydrodynamical 
simulations have convincingly shown that the IGM 
does fragment into structures at early times in CDM 
cosmogonies (e.g. Cen et al.l994; Zhang et al.l995; 
Hernquist et al.l996), the same simulations are much 
less able to predict the efficiency with which the first 
gravitationally collapsed objects lit up the universe 
at the end of the dark age. 

The scenario that has received the most theoret- 
ical studies is one where hydrogen is photoionized 
by the UV radiation emitted either by quasars or by 
stars with masses > 10 M© (e.g. Shapiro & Giroux 
1987; Haiman & Loeb 1998; Madau et al.l999; Chiu 
& Ostriker 2000; Ciardi et al.2000), rather than ion- 
ized by collisions with electrons heated up by, e.g. 



supernova-driven winds from early pregalactic ob- 
jects. In the former case a high degree of ionization 
requires about 13.6 x (1 + t/frec)eV per hydrogen 
atom, where free is the volume-averaged hydrogen 
recombination timescale, t/free being much greater 
than unity already at z ~ 10 according to numer- 
ical simulations. Collisional ionization to a neutral 
fraction of only a few parts in 10^ requires a compa- 
rable energy input, i. e. an IGM temperature close 
to 10^ K or about 25 eV per atom. 

Massive stars will deposit both radiative and me- 
chanical energy into the interstellar medium of pro- 
togalaxies. A complex network of 'feedback' mecha- 
nisms is likely at work in these systems, as the gas in 
shallow potential is more easily blown away (Dekel 
& Silk 1986; Tegmark et al.l993; Mac Low & Fer- 
rara 1999; Mori et al. 2002) thereby quenching star 
formation. Furthermore, as the blastwaves produced 
by supernova explosions reheat the surrounding in- 
tergalactic gas and enrich it with newly formed heavy 
elements (see below), they can inhibit the formation 
of surrounding low-mass galaxies due to 'baryonic 
stripping' (Scannapieco et al.2002). It is therefore 
difficult to establish whether an early input of me- 
chanical energy will actually play a major role in 
determining the thermal and ionization state of the 
IGM on large scales. What can be easily shown is 
that, during the evolution of a a 'typical' stellar pop- 
ulation, more energy is lost in ultraviolet radiation 
than in mechanical form. This is because in nuclear 
burning from zero to solar metallicity {Zq = 0.02), 
the energy radiated per baryon is 0.02 x 0.007 x muc^ , 
with about one third of it going into H-ionizing pho- 
tons. The same massive stars that dominate the 
UV light also explode as supernovae (SNe), return- 
ing most of the metals to the interstellar medium 
and injecting about 10^^ ergs per event in kinetic en- 
ergy. For a Salpeter initial mass function (IMF), 
one has about one SN every 150 Mq of baryons that 
forms stars. The mass fraction in mechanical energy 
is then approximately 4 x 10~^, ten times lower than 
the fraction released in photons above 1 ryd. 

The relative importance of photoionization ver- 
sus shock ionization will depend, however, on the 
efficiency with which radiation and mechanical en- 
ergy actually escape into the IGM. Consider, for ex- 
ample, the case of an early generation of subgalac- 
tic systems collapsing at redshift 9 from 2-cr fluc- 
tuations. At these epochs their dark matter halos 
would have virial radii rvir = 0.75h~^ Kpc and cir- 
cular velocities Vc(^vir) = 25 kms~^, corresponding 
in top-hat spherical collapse to a virial tempera- 
ture Tvir = 0.5jj.mpV^ /k ~ 10^"^ K and halo mass 
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M = O.IV^ /GH ^ lO^h-^ M0.2 Halos in this 
mass range are characterized by very short dynam- 
ical timescales (and even shorter gas coohng times 
due to atomic hydrogen) and may therefore form 
stars in a rapid but intense burst before SN 'feed- 
back' quenches further star formation. For a star for- 
mation efficiency of / = 0.1, Q^h'^ = 0.02^, h = 0.5, 
and a Salpeter IMF, the explosive output of 10, 000 
SNe will inject an energy ^ 10^^ ergs. This is 
roughly a hundred times higher than the gas binding 
energy: a significant fraction of the halo gas will then 
be lifted out of the potential well ('blow-away') and 
shock the intergalactic medium (Madau et al.2001). 
If the explosion occurs at cosmic time t = 4 x 10^ yr 
(corresponding in the adopted cosmology to 2: = 9), 
at time At = OAt after the event it is a good ap- 
proximation to treat the cosmological blast wave as 
adiabatic, with proper radius given by the standard 
Sedov-Taylor self-similar solution, 

i?.«(i^^)''%2/^At2/^«17kpc. (1) 

Here rj ^ 0.3 is the fraction of the available SN 
energy that is converted into kinetic energy of the 
blown-away material (Mori et al.2002). At this in- 
stant the shock velocity relative to the Hubble fiow 
is 

Vs ^ 2Rs/bAt ^ 40 kms"\ (2) 

lower than the escape velocity from the halo cen- 
ter. The gas temperature just behind the shock front 
is Tg = 3jj.mpVg /16k > 10^*^ K, enough to ionize 
all incoming intergalactic hydrogen. At these red- 
shifts, it is the onset of Compton cooling off cos- 
mic microwave background photons that ends the 
adiabatic stage of blast wave propagation; the shell 
then enters the snowplough phase and is finally con- 
fined by the IGM pressure. According to the Press- 
Schechter formalism, the comoving abundance of col- 
lapsed dark halos with mass M = 10^ Mq at 
z = 9 is dn/dlnM ~ 80/i^Mpc~^, corresponding 
to a mean proper distance between neighboring ha- 
los of ^ lb h~^kpc. With the assumed star forma- 
tion efficiency, only a small fraction, about 4 percent, 
of the total stellar mass inferred today (Fukugita et 
al.l998) would actually form at these early epochs. 
Still, our simple analysis shows that the blast waves 
from such a population of pregalactic objects could 
drive vast portions of the IGM to a significantly 

^This assumes an Einstein-de Sitter universe with Hubble 
constant Hq = 100 h kms"-*^ Mpc"-*^. 

^Here 0,^, is the baryon density parameter, and ffli) is the 
fraction of halo mass converted into stars. 



higher adiabat, Tigm ~ 10^ K, than expected from 
photoionization, so as to 'choke off' the collapse of 
further M < 10^ Mq systems by raising the cos- 
mological Jeans mass. In this sense the process may 
be self-regulating. 

The thermal history of expanding intergalactic 
primordial gas at the mean density is plotted in Fig- 
ure 1 as a function of redshift for a number of illus- 
trative cases. The code we have used includes the 
relevant cooling and heating processes and follows 
the non-equilibrium evolution of hydrogen and he- 
lium ionic species in a cosmological context. 




Fig. 1: Thermal history of intergalactic gas at the mean 
density in an Einstein-de Sitter universe with l^fo/i^ = 0.02 
and h = 0.5. Short dash-dotted line: temperature evolution 
when the only heating source is a constant ultraviolet (CUV) 
background of intensity 10~^^ ergs cm~^ s"-*^ Hz"-*^ sr~^ at 1 
Ryd and power-law spectrum with energy slope a = 1. Long 
dash-dotted line: same for the time-dependent quasar ion- 
izing background as computed by Haardt &; Madau (1996; 
HM). Short-dashed line: heating due to a CUV background 
but with an initial temperature of 4 x 10^ K at 2; = 9 as ex- 
pected from an early era of pregalactic outflows. Long-dashed 
line: same but for a HM background. Solid line: heating due 
to a CUV background but with an initial temperature of 10^ K 
at 2; = 9. Dotted line: same but for a HM background. 

The gas is allowed to interact with the CMB 
through Compton cooling and either with a 
time-dependent QSO ionizing background as com- 
puted by Haardt & Madau (1996) or with a 
time-independent metagalactic flux of intensity 
10~^^ ergscm"^ s~^ Hz~^ sr~^ at 1 Ryd (and power- 
law spectrum with energy slope a = 1). The tem- 
perature of the medium at z = 9 - where we start 
our integration - has been either computed self- 
consistently from photoheating or fixed to be in the 
range 10^"^ — 10^ K expected from SN-driven bubbles 
with significant filling factors (see below). The var- 
ious curves show that the temperature of the IGM 
at z = 3 — 4 will retain little memory of an early 
era of pregalactic outfiows and preheating, and be 
consistent with that expected from photoionization. 
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3. FEEDBACK AND PREGALACTIC 
ENRICHMENT 

Understanding the origin of the chemical ele- 
ments, following the increase in their abundances 
with cosmic time, and uncovering the processes re- 
sponsible for distributing the products of stellar nu- 
cleosynthesis over very large distances are all key 
aspects of the evolution of gaseous matter in the 
universe. One of the major discoveries with Keck 
concerning the IGM has been the identification 
of metal absorption lines associated with many of 
the Lya forest systems. The detection of measur- 
able C IV and Si iv absorption lines in clouds with 
H I column densities as low as 10^^"^ cm~^ implies 
a minimum universal metallicity relative to solar in 
the range [-3.2] to [-2.5] at z - 3 (Songaila 1997). 
There is no indication in the data of a turnover 
in the C iv column density distribution down to 
A^civ < 10^^ '^ cm-2 (TVhi < lO^'^-^ cm-^, Ellison 
et al.2000). 

From a theoretical perspective, it is unclear 
whether the existence of heavy elements in the IGM 
at z = 3 — 3.5 points to an early {z > 6) enrichment 
epoch by low-mass subgalactic systems (Madau et 
al.2001), or is rather due to late pollution by the pop- 
ulation of star-forming galaxies known to be already 
in place at z = 3. The Press-Schechter theory for the 
evolving mass function of dark matter halos predicts 
a power-law dependence, dN/d\nm ex TT^^^^ff ~^)/^, 
where ngff is the effective slope of the CDM power 
spectrum, ngff ~ —2.5 on subgalactic scales. As 
hot, metal-enriched gas from SN-driven winds es- 
capes its host halo, shocks the IGM, and eventually 
forms a blast wave, it sweeps a region of intergalac- 
tic space the volume of which increases with the 3/5 
power of the injected energy Eq (in the adiabatic 
Sedov-Taylor phase). The total fractional volume 
or porosity, Q, filled by these 'metal bubbles' per 
unit explosive energy density Eq dN / dlnm is then 
Q (X E^^^dN/d\nm oc (dTV/dlnm)^/^ ex m-^^/^^ 
Within this simple scenario it is the star-forming 
objects with the smallest masses which will arguably 
be the most efficient pollutant of the IGM on large 
scales. Metal-enriched material from SN ejecta may 
be far more easily accelerated to escape velocities in 
the shallow potential wells of subgalactic systems at 
z > 6. Late enrichment may also encounter problems 
in explaining the kinematic quiescence of C iv lines; 
the observed small scale properties of the IGM at 
z ^ 3 appear consistent with C iv absorbers be- 
ing the result of ancient pregalactic outflows (Ranch 
et al.2001). According to numerical hydrodynamics 
simulations of structure formation in the IGM, the 



metals associated with logA^Hi ^ 14 filaments fill a 
fraction > 5% of intergalactic space, and are there- 
fore far away from the high overdensity peaks where 
galaxies form, gas cools, and star formation takes 
place. Their chemical enrichment may then reflect 
more uniform (i.e. 'early') rather than in-situ (i.e. 
'late') metal pollution. The case for pregalactic en- 
richment may have been recently strenghtened by 
the observation of an invariant C iv column density 
distribution throughout the redshift range 2 < z < 5 
(Songaila 2001). 
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Fig. 2: Snapshots of the logarithmic number density of the 
gas at five different elapsed times for our Case 1 simulation 
run. The three panels in each row show the spatial density 
distribution in the X — Y plane on the nested grids. The three 
columns in each figure depict the time evolution from about 
6 Myr to up to 180 Myr. Along a given row, the leftmost 
panel refers to grid L5 (linear size 48 kpc), the central one to 
grid L3 (linear size 12 kpc), and the rightmost panel refers to 
the grid LI (linear size 3 kpc). The density range is —5 < 
log (n/cm-3) < 1. The halo gas is assumed to be initially in 
hydrostatic equilibrium and non self-gravitating. (From Mori 
et al.2002.) 

To simulate the process of blow-away we have de- 
veloped a 3D Eulerian code that solves the hydro- 
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dynamic equations for a perfect fluid in Cartesian 
geometry. To deal with very different length scales 
in our simulation we have adopted a 'nested grid 
method' with six levels of fixed Cartesian grids. The 
grids are connected by the transfer of conserved vari- 
ables, and are centered within each other, with the 
finest covering the whole galaxy halo. Since the cell 
number is the same (128 x 128 x 128) for every level 
Ln {n = 1,2,..., 6), the minimum resolved scale is 
about 22 pc and the size of the coarsest grid is 96 
kpc. Thus, the scheme has a wide dynamic range in 
the space dimension. 

The results of a numerical simulation (run on 
massive parallel supercomputers at the Center for 
Computational Physics, Tsukuba University) are de- 
picted in Figure 2. In this run about 10,000 SNe 
explode in a high redshift protohalo (correspond- 
ing to a star formation efficiency of about 10% for 
a Salpeter IMF). Our algorithm for simulating SN 
feedback improves upon previous treatments in sev- 
eral ways. OB associations are distributed as a func- 
tion of gas density according to a Schmidt-type law 
((X p^) using a Monte-Carlo procedure. After a main 
sequence lifetime, all stars more massive than 8 Mq 
explode instantaneously injecting an energy of 10^^ 
ergs, and their outer layers are blown out leaving a 
compact remnant of 1.4 Mq. Therefore SNe inject 
energy (assumed in pure thermal form) and mass 
into the interstellar medium: these are supplied to a 
sphere corresponding to the radius of a SN remnant 
in a uniform ambient medium of density n and in 
the adiabatic Sedov-Taylor phase; the expansion ve- 
locity is also self-consistently calculated from such 
solution. The gas then starts cooling immediately 
according to the adopted gas cooling function. While 
the time-step At is controlled by the Courant con- 
dition, if there are more than two SN events in a OB 
association we decrease At until only one explosion 
per association occurs during the time-step. 

Snapshots of the gas density distribution in a pla- 
nar slice of the nested grids are shown for an ex- 
tended stellar distribution case {a = 1). After a few 
Myr from the beginning of the simulation the most 
massive stars explode as SNe and produce expand- 
ing hot bubbles surrounded by a cooling dense (n ^ 
1 cm~^) shell. As the evolution continues, a coher- 
ent and increasingly spherical shell expanding into 
the IGM is eventually formed. The shell contains a 
large fraction of the halo gas that has been swept-up 
during the evolution. The final bottom row of the 
simulation figure shows the final stages of the evolu- 
tion: the shell is now nearly spherically symmetric, 
its interior being filled with warm (T < 10^ K) gas 



at a very low density n < 10~^ cm~^, i.e. slightly 
below the mean value for the IGM. At the end of 
the simulations, the shell is still sweeping out IGM 
material; its radius and velocity are 21 kpc and 26 
km sec~^ at t = 250 Myr. Using momentum conser- 
vation one can estimate the final radius of the shell 
to be close to 25 kpc. This is about 17 times the 
virial radius of the halo! 

It is clear then that SN-driven pregalactic out- 
fiows may be an efficient mechanism for spreading 
metals around. The collective explosive output of 
about ten thousands SNe per M > 10^ h'^ M© 
halo at these early epochs could then pollute the en- 
tire intergalactic space to a mean metallicity (Z) = 
^z/^b ^ 0.003 (comparable to the levels observed 
in the Lya forest at 2; ~ 3) without much perturb- 
ing the IGM hydrodynamically, i.e. producing large 
variations of the baryons relative to the dark mat- 
ter. The significance of these results goes well be- 
yond early metal enrichment. This is because, since 
the cooling time of collisionally ionized high density 
gas in small halos at high redshifts is much shorter 
than the then Hubble time, virtually all baryons are 
predicted to sink to the centers of these halos in 
the absence of any countervailing effect (White & 
Rees 1978). Efficient feedback is then necessary in 
hierarchical clustering scenarios to avoid this 'cool- 
ing catastrophe', i.e. to prevent too many baryons 
from turning into stars as soon as the first levels of 
the hierarchy collapse. The required reduction of the 
stellar birthrate in halos with low circular velocities 
may naturally result from the heating and expulsion 
of material due to OB stellar winds and repeated SN 
explosions from the first burst of star formation. 

4. PROBING THE EPOCH OF FIRST LIGHT AT 

21-CM 

Prior to the epoch of full reionization, the inter- 
galactic medium and gravitationally collapsed sys- 
tems may be detectable in emission or absorption 
against the CMB at the frequency corresponding to 
the redshifted 21-cm line (associated with the spin- 
fiip transition from the triplet to the singlet state of 
neutral hydrogen.) In general, 21-cm spectral fea- 
tures are expected to display angular structure as 
well as structure in redshift space due to inhomo- 
geneities in the gas density field, hydrogen ionized 
fraction, and spin temperature. Several different sig- 
natures have been investigated in the recent litera- 
ture: (a) the ffuctuations in the redshifted 21-cm 
emission induced by the gas density inhomogeneities 
that develop at early times in CDM-dominated cos- 
mologies (Madau et al.l997; Tozzi et al.2000) and 
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by virialized "minihalos" with Tvir < 10^ K (Iliev et 
al.2002); (b) the sharp absorption feature in the ra- 
dio sky due to the rapid rise of the Lya continuum 
background that marks the birth of the first UV 
sources in the universe (Shaver et al.l999); (c) the 
21-cm narrow hues generated in absorption against 
very high-redshift radio sources by the neutral IGM 
(Carihi et al.2002) and by intervening minihalos and 
protogalactic disks (Furlanetto & Loeb 2002). 

A quick summary of the physics of 21-cm radi- 
ation will illustrate the basic ideas and unresolved 
issues behind these studies. The emission or absorp- 
tion of 21-cm photons from a neutral IGM is gov- 
erned by the hydrogen spin temperature Ts defined 
by ni/no = 3 exp(— T^/T^), where no and ni are 
the singlet and triplet n = 1 hyperfine levels and 
ksT^ = 5.9 X 10~^eV is the energy of the 21-cm 
transition. In the presence of only the CMB radi- 
ation with TcMB = 2.73 (1 + z) K, the spin states 
will reach thermal equilibrium with the CMB on a 
timescale of T*/(Tcmb^io) ^ 3 x 10^ (1 + z)-^yi 
(Aio = 2.9 X 10~^^ s~^ is the spontaneous decay rate 
of the hyperfine transition of atomic hydrogen), and 
intergalactic H i will produce neither an absorption 
nor an emission signature. A mechanism is required 
that decouples Ts and Tcmb, e.g. by coupling the 
spin temperature instead to the kinetic temperature 
Tk of the gas itself. Two mechanisms are avail- 
able, collisions between hydrogen atoms (Purcell & 
Field 1956) and scattering by Lya photons (Field 
1958). The collision-induced coupling between the 
spin and kinetic temperatures is dominated by the 
spin-exchange process between the colliding hydro- 
gen atoms. The rate, however, is too small for re- 
alistic IGM densities at the redshifts of interest, al- 
though collisions may be important in dense regions 
with dp/ p ^ 30[(1 + 2;)/10]~^, like virialized miniha- 
los. 

In the low density IGM, the dominant mechanism 
is the scattering of continuum UV photons redshifted 
by the Hubble expansion into local Lya photons. 
The many scatterings mix the hyperfine levels of 
neutral hydrogen in its ground state via interme- 
diate transitions to the 2p state, the Wouthuy sen- 
Field process. As the neutral IGM is highly opaque 
to resonant scattering, the shape of the continuum 
radiation spectrum near Lya will follow a Boltz- 
mann distribution with a temperature given by the 
kinetic temperature of the IGM (Field 1959). In this 
case the spin temperature of neutral hydrogen is a 
weighted mean between the matter and CMB tem- 
peratures. There exists then a critical value of the 
background flux of Lya photons which, if greatly ex- 



ceeded, would drive the spin temperature away from 

^CMB- 

While the microphysics is well understood, our 
understanding of the astrophysics of 21-cm tomog- 
raphy is still poor. In Tozzi et al.(2000) we used 
N-body cosmological simulations and, assuming a 
fully neutral medium with Ts ^ Tcmb, showed that 
prior to reionization the same network of sheets and 
filaments (the 'cosmic web') that gives rise to the 
Lya forest at 2: ~ 3 should lead to fluctuations in 
the 21-cm brightness temperature at higher redshifts 
(Figure 3). At 150 MHz {z = 8.5), for observations 
with a bandwidth of 1 MHz, the root mean square 
fluctuations should be ~ 10 mK at 1', decreasing 
with scale. Because of the smoothness of the CMB 
sky, fluctuations in the 21-cm radiation will dom- 
inate the CMB fluctuations by about 2 orders of 
magnitude on arcmin scales. The search at 21-cm 
for the epoch of first light has become one of the 
main science drivers of the LOw Frequency ARray 
{LOFAR] see http://www.astron.nl/lofar/science/). 
While remaining an extremely challenging project 
due to foreground contamination from extragalactic 
radio sources (Di Matteo et al.2002), the detection 
and imaging of these small-scale structures with LO- 
FAR is a tantalizing possibility within range of the 
thermal noise of the array. 
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Fig. 3: Left: Radio map of redshifted 21-cm emission against 
the CMB in a 'tilted' CDM (tCDM) cosmology at z = 8.5. 
Here a collisionless A^-body simulation with 64^ particles has 
been performed with Hydra (Couchman et al.l995). The sim- 
ulation box size is 20h~^ comoving Mpc, corresponding to 17 
(11) arcmin in tCDM (OCDM). The baryons are assumed to 
trace the dark matter distribution without any biasing, the 
spin temperature to be much greater than the temperature of 
the CMB everywhere, and the gas to be fully neutral. The 
point spread function of the synthesized telescope beam is a 
spherical top-hat with a width of 2 arcmin. The frequency 
window is 1 MHz around a central frequency of 150 MHz. 
The contour levels outline regions with signal greater than 
4/iJy per beam. Right: Same for a open cosmology (OCDM). 
Since the growth of density fluctuations ceases early on in an 
open universe (and the power spectrum is normalized to the 
abundance of clusters today), the signal at a given angular 
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size is much larger in OCDM than in tCDM at these early 
epochs. (From Tozzi et al.2000.) 

On the theoretical side, there are several effects that 
need to be examined further. As mentioned above, 
it is the presence of a sufficient flux of Lya photons 
which renders the neutral IGM 'visible'. Without 
heating sources, the adiabatic expansion of the uni- 
verse will lower the kinetic temperature of the gas 
well below that of the CMB, and the IGM will be de- 
tectable through its absorption. If there are sources 
of radiation that preheat the IGM, it may be pos- 
sible to detect it in emission instead. The ener- 
getic demand for heating the IGM above the CMB 
temperature is meager, only ~ 0.004 eV per parti- 
cle at z ~ 10. Consequently, even relatively ineffi- 
cient heating mechanisms may be important warm- 
ing sources well before the universe was actually 
reionized. Perhaps more importantly, prior to full 
reionization the IGM will be a mixture of neutral, 
partially ionized, and fully ionized structures: low- 
density regions will be fully ionized first, followed 
by regions with higher and higher densities. Radio 
maps at 21 {l-\-z) cm will show a patchwork of emis- 
sion/absorption signals from H i zones modulated by 
H II regions where no signal is detectable against the 
CMB. It is the early generation of stars likely re- 
sponsible for reionization which will also generate a 
background radiation field of photons with energies 
between 10.2-13.6 eV to which the IGM is transpar- 
ent. As each of these photons gets redshifted, it ulti- 
mately reaches the Lya transition energy of 10.2 eV, 
scatters resonantly off neutral hydrogen, and mixes 
the hyperfine levels. The observability of the pre- 
reionization IGM depends in this case on the UV 
spectrum of the firsts stars, i.e. on the number of 
Lya continuum photons emitted per H-ionizing pho- 
ton. 

The integrated UV spectrum of a stellar popu- 
lation is characterized by a strong break at the Ly- 
man edge whose size depends on age, initial mass 
function, and star formation history. In Figure 4 we 
have therefore normalized the Lya rate to the rate 
of emission of hydrogen-ionizing photons; after inte- 
grating over the stellar age, r, a starburst is char- 
acterized by a ratio TVa/A^ion ^ 22 (r/0.1 Gyr)^-^^ 
(Scalo IMF), about three times higher than the con- 
stant SFR case (this is because of the shorter lifetime 
of the massive stars that produce 1 ryd radiation). 
It is possible that the startburst mode may be more 
relevant for the galaxies responsible for reionization, 
as stellar winds and supernovae can easily expell the 
gas out of the shallow potential wells of low mass ha- 



los, after the first burst of star formation occurs. In 
this case the escape fraction might be low until the 
gas is expelled, and as Figure 4 indicates this would 
contribute to an even stronger Lya continuum emis- 
sion per ionizing flux. 




log age (years) log age (years) 

Fig. 4: Left: Continuum Lyo; photon emission rate Na 
(s"-*^) versus age (in years) for a stellar population with a 
Salpeter IMF (dashed lines) and a Scalo IMF (solid lines). 
Upper curves assume an instantaneous burst of star forma- 
tion, lower curves a constant star formation rate. The rate of 
production of Lyo; photons has been normalized to the rate 
of emission of hydrogen-ionizing photons, iVion- The popula- 
tion synthesis values have been computed for low metallicities 
(Z = 0.02 Zq), and are based on an update of Bruzual & 
Chariot's (1993) libraries. Right: Same after integration over 
the age of the stellar population. 



5. PHOTON CONSUMPTION IN MINIHALOS 
DURING COSMOLOGICAL REIONIZATION 

The process of reionization has recently received 
much theoretical attention, and has been studied 
by several authors using semi-analytic models (e.g. 
Meiksin & Madau 1993; Shapiro et al.l994; Haiman 
& Loeb 1997, 1998; Benson et al.2001) and three- 
dimensional numerical simulations (Gnedin 2000). 
In general, these works aim to follow the time evo- 
lution of the filling factor of ionized (H ii ) regions, 
based on some input prescriptions for the emissivity 
and spectra of the ionizing sources. 

More recent works have focused on the increased 
rate of recombinations in a clumpy medium relative 
to a homogeneous one, i.e. when (p^) > (p)^ (Cia- 
rdi et al. 2000; Benson et al.2001; Chiu & Ostriker 
2000; Gnedin 2000; Madau et al.l999). These stud- 
ies have left significant uncertainties on the details 
of how reionization proceeds in an inhomogeneous 
medium. Since the ionizing sources are embedded 
in dense regions, one might expect that these dense 
regions are ionized first, before the radiation escapes 
to ionize the low-density IGM. Alternatively, most of 
the radiation might escape from the local, dense re- 
gions along low column density lines of sight. In this 



8 



MADAU 



case, the under dense 'voids' are ionized first, with 
the ionization of the denser filaments and halos lag- 
ging behind (Miralda-Escude et al. 2000). 

At the earliest epochs of structure formation in 
CDM cosmologies the smallest nonlinear objects are 
the numerous small halos that condense with virial 
temperatures between the cosmological Jean tem- 
perature and ~ 10^ K. Such "minihalos" are not 
yet fully resolved nor is the process of their photoe- 
vaporation captured in large-scale three-dimensional 
cosmological simulations.^ In Haiman et al.(2001) 
we have used semi-analytic methods combined with 
three-dimensional numerical simulations of individ- 
ual photoevaporating minihalos to (1) quantify the 
importance of high-redshift minihalos as sinks of ion- 
izing photons; and (2) assess whether by extrapolat- 
ing to early times the known population of galaxies 
and quasars a sufficient number of UV photons are 
produced for hydrogen reionization. If reionization 
occurs at sufficiently high redshifts {zr ^ 20), the 
intergalactic medium is heated to ~ lO^K and most 
minihalos never form. 




Fig. 5: Left: The evolution of the density profile of a mini- 
halo with mass 10^ M©, illuminated by a UV background fiux 
at z = 10. The hydrogen photoionization rate is assumed to 
be uniform on the grid, i.e. we do not solve radiative trans- 
fer. Photoevaporation occurs very rapidly, with most recom- 
binations occurring within a sound-crossing time. The figure 
reveals that by 2; = 9.5 the density contrast is everywhere re- 
duced below 6 ^ 10. Right: Total number of recombinations 
inside minihalos, following sudden reionization at redshift Zr- 
Recombinations are shown per background hydrogen atom. 
In the earliest stages of reionization, minihalos dominate the 
average gas clumping in the universe, and can therefore domi- 
nate the total recombination rate. (From Haiman et al.2001.) 

On the other hand, if Zr ^ 20 as it now appears 

^Photoionization by either a nearby external UV source, 
or by the smooth cosmic UV background radiation after the 
reionization epoch will heat the gas inside a minihalo to a 
temperature ~ 10^ K. By definition, the gas inside a minihalo 
is then no longer bound, leading to the photoevaporation of 
baryons out of their host halos (Shapiro, Raga & Mellema 
1998; Barkana & Loeb 1999). 



more likely, then a significant fraction 10%) of all 
baryons have already collapsed into minihalos, and 
are subsequently removed from the halos by photo- 
evaporation as the ionizing background flux builds 
up (Figure 5). This process requires a significant 
budget of ionizing photons, about 10-20 of them per 
background hydrogen atom. This exceeds the pro- 
duction by a straightforward extrapolation back in 
time of known quasar and galaxy populations by a 
factor of up to ~ 10 and ~ 3, respectively (Haiman 
et al.2001). 

6. THE EARLIEST LUMINOUS SOURCES AND 
THE WING OF THE GUNN-PETERSON 
TROUGH 

Prior to complete reionization at redshift z^, 
sources of ultraviolet radiation will be seen behind 
a large column of intervening gas that is still neu- 
tral. In this case, because of scattering off the line- 
of-sight due to the diffuse neutral IGM, the spectrum 
of a source at Zem > should show the red damping 
wing of the Gunn-Peterson (1965) absorption trough 
at wavelengths longer than the local Lya resonance, 
Aobs > Ac(l + 2;em), whcrc \a, = c/va = 1216 A. 
At ^ 6, this characteristic feature extends for 
more than 1500 kms~^ to the red of the resonance, 
and may significantly suppress the hya emission line 
in the spectra of the first generation of objects in 
the universe. Measuring the shape of the absorption 
profile of the damping wing could provide a deter- 
mination of the density of the neutral IGM near the 
source (Miralda-Escude 1998). 

A number of authors (Madau & Rees 2000; Gen 
& Haiman 2000; Haiman 2002) have recently focused 
on the width of the red damping wing - related to 
the expected strength of the Lya emission line - in 
the spectra of very distant QSOs and star-forming 
galaxies as a flag of the observation of the IGM 
before reionization. We have assessed, in particu- 
lar, the impact of the photoionized, Mpc-size regions 
which will surround individual luminous sources of 
Lyman continuum radiation on the transmission of 
photons redward of the Lya resonance, and shown 
that the damping wing of the Gunn-Peterson trough 
may nearly completely disappear because of the lack 
of neutral hydrogen in the vicinity of a bright QSO. 
The effect of this local photoionization is to greatly 
reduce the scattering opacity between the redshift of 
the quasar and the boundary of its H ii region (Fig- 
ure 6). 

The detection of a strong Lya emission line in 
the spectra of bright QSOs shining for > 10^ yr can- 
not then be used, by itself, as a constraint on the 
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reionization epoch. The first signs of an object radi- 
ating prior to the transition from a neutral to an 
ionized universe may be best searched for in the 
spectra of luminous sources with a small escape frac- 
tion of Lyman continuum photons into the IGM, or 
sources with a short duty cycle (like, e.g., gamma-ray 
bursts). 




9300 9600 9900 9300 9600 9900 
observed wavelength (A) 

Fig. 6: The Keck/LRIS spectrum of the faint Zem = 5.5 
quasar RD J0301117+002025 (Stern et al.2000), redshifted to 
Zem = 7. The dotted vertical Hues indicate the expected wave- 
length of the Lya resonance. Fluxes have arbitrary normaliza- 
tions. The dashed curves show the transmission exp(— tgp) 
(where tgp is the Gunn-Peterson optical depth to resonant 
scattering) through a uniform IGM, assuming the QSO is be- 
ing observed prior to the reionization epoch at Zr = Q. Four 
different cases are shown, as the rate of emission of Lyman 
continuum photons which ionize the IGM in the vicinity of 
the source is varied from Ni = to = 10^^, lO^'^, and 10^^ 



The calculations assume a quasar lifetime of ts = 10^ 



a power-law spectrum of the form fu oc with a = 0.5 

near the hydrogen Lyman edge, a recombination timescale of 
tree = 1.3 Gyr, an Einstein-de Sitter universe with h = 0.5, 
and a baryon density parameter Qsh^ = 0.02. (From Madau 
& Rees 2000.) 
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